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' Abstract. This work focuses on the properties of dusty tori in active galactic 

\ nuclei (AGN) derived from the comparison of SDSS type 1 quasars with mid- 

Infrared (MIR) counterparts and a new, detailed torus model. The infrared data 
D ' were taken by the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. 

^0 ! Basic model parameters are constraint, such as the density law of the graphite 

and silicate grains, the torus size and its opening angle. A whole variety of 
optical depths is supported. The favoured models are those with decreasing 
density with distance from the centre, while there is no clear tendency as to 
the covering factor, i.e. small, medium and large covering factors are almost 
equally distributed. Based on the models that better describe the observed 
fS| , SEDs, properties such as the accretion luminosity, the mass of dust, the inner 

' to outer radius ratio and the hydrogen column density are computed. The 

0^ ' properties of the tori, as derived fitting the observed SEDs, are independent of 

, the redshift, once observational biases are taken into account. 
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1. Torus Model 



o 

c/3 . A detailed description of the torus model is given in iFritz et all (|2nn(J l: here 

^ ' we briefly present the main characteristics. The adopted torus geometry is a 

flared disk, i.e. a sphere with the polar cones removed. Its size is defined by its 
outer radius, Rout, and its opening angle. The dust components that dominate 
both the absorption and the emission of radiation are graphite and silicate. The 
inner radius, Rin, depends both on the sublimation temperature of the dust 
grains (1500 and 1000 K, for graphite and silicate, respectively) and on the 
accretion luminosity. Since silicate grains have a lower sublimation temperature 
than graphite grains, the innermost regions of the torus only cons ist of graphite. 
The a dopted absorption and scattering coefficients are those bv lLaor &: DrainI 
for d ust grains of difFeren t dimensions, weighted with the standard MRN 
distribution ( Mathis et al.lll97'^ ) . The gas density within the torus is modeled in 



such a way to allow a gradient along both the radial and the angular coordinates. 

The central source is assumed to be point-like and its emission isotropic. Its 
spectral energy distribution is defined by means of a composition of power laws 
with different values for the spectra index in the UV, optical and IR. The A- 
iteration method is used to solve the radiative transfer equation. A geometrical 
grid is defined along the three spatial coordinates, and the main physical quan- 
tities (dust density and temperatures, electromagnetic emission, optical depth, 
etc.) are computed with respect to the center of the volume elements defined by 
the grid. The global SED is computed at difFerent angles of the line-of-sight with 
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respect to the torus' equatorial plane and includes three contributions: emission 
from the AGN, thermal emission and scattering emission by dust in each volume 
element. An example of an emitted spectrum is given in Fig. ^ 
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Figure 1. Emission spectra as a function of wavelength for 10 different lines- 
of-sight inclinations from 5" = 0° (edge-on; lower curve) to ^' = 90° (face-on, 
upper curve) at regular steps of 10°, for a geometrical configuration with 
r(9.7) = 10.0, Rout/Rin — 300, torus opening angle 6 — 140° and density 
law p (r, 0) oc 7-~^e^^l^°'''^^^l, where r is the distance from the centre and 9 the 
angle from the equatorial plane. 



2. Data 

The MIR data used here are taken from the SWIRE ENl and EN2 and Lockman 
fields and the SWIRE catalogues p rocessed by the SWIRE collaboration. De - 
tails about the data can be found in lLonsdale et ah |(|200l),|S urace et ahl (|200l) 



and IShupe et all (|2006|). The optical data were taken from the Sloan Digital 
Sky Survey (SDSS) Data Release 4 (DR4), that covers the entire Lockman and 
SWIRE EN2 fields and a part of the SWIRE ENl field. 2MASS data points 
were also use d, whenever available as well as data from the Deep 2MASS Survey 
in Lockman ( Beichman et al.ll2003r ). 



A total of 279 spectroscopically confirmed quasars, all detected in all IRAC 
bands and MIPS 24 micron, lie within the SWIRE fields covered by the SDSS 
DR4 spectroscopic release, with redshifts spanning from 0.214 to 5.215. Their i- 
band magnitudes reach 19.1 for objects with redshifts typically less than 2.3 and 
go up to a magnitude deeper for higher redshifts (Richard s et ah 200^. For a 
detailed analysis of the properties of the sample in ENl see lHatziminaoglou et all 
(2003). 
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3. First Results and Implications 



Each observed SED is compared to a total of 720 models, with 10 lines of sight 
each (from an equatorial to a pole-on view). Examples of best fits are shown in 
Fig. 121 Some of the results of the SED fitting are summarised below. 
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Figure 2. Examples of SED fits for three AGN of our sample. The fluxes 
are given in Jy, the wavelengths in A. 



For ~70% of the objects the density, p{r,9) oc r^e"'^''^"''^^^!, decreases radially 
from the centre {(3 < 0) while for another ~25% it remains constant (/3 = 0). 
Half of the objects have a density distribution that does not vary with the 
angle, 6, from the equator. Trying not to exclude a priori any of the model 
characteristics, we allow for low optical depth tori. Some 65% of the objects were 
better described by models with rg.y < 1.0. We therefore need to recalculate 
the covering factor and consider that we switch from a type 1 to a type 2 object 
(i.e. that we no longer see the central source through the torus) when tq.s > 1. 

Inner radii are typically of order of a few pc, with extreme cases (very 
luminous objects) reaching ~ 10 pc, as seen on the left panel of Fig. |3J The 
mass of dust, MDust) is computed summing the individual sample elements of 
the best-fit model. Typical values range from 10^ to 10*^ Mq. Note that Moust 
does not refer to the total mass of the torus, that can be obtained by adding the 
mass of gas, typically ~100 times larger than that of the dust. Many quantities 
(Rin, Rout! MDust show a clearly increasing trend with redshift, z. This could 
be interpreted as an indication for larger and/or more massive tori at higher 
z; however it is simply the manifestation of the Scott effect, i.e. the fact that 
in order for distant objects to make it into a flux-limited sample they need to 
be intrinsically brighter. What happens in reality is that the inner radius, Rin, 
directly scales with the accretion luminosity. Lace (expressed in units of 10"^^ 
erg/sec), via the relation bv iBarvainid ( 19871 1: Rm — 1.3 x y^Lacc) x '^lioo [P^]) 
where the temperature, T, is given in units of 1500K. In a flux-limited sample 
Lace; which is directly related to the bolometric luminosity, increases with z and 
therefore so does Rin, as well as Rout and MDust that directly depend on Rin, 
as already mentioned. We therefore deduce that sources of the same accretion 
luminosity can have similar properties, independently on the redshift. 

The eg'watona/ hydrogen column density, Nh, distribution, derived from the 
models, is shown in the right panel of Fig. |21 The values are representative 
of type 2 objects, with some in the range of Compton thick AGN (Nh ~ 10^^ 
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Figure 3. Left panel: histogram of the inner torus radii, Rjn; Right panel: 
distribution of the hydrogen column density, Nh, along the equator. 

cm~^). If the distinction between type 1 and type 2 objects is simply a matter of 
orientation, as suggested from the Unified Scheme, this would be the expected 
distribution of equatorial Nh- 

4. Conclusions 

This work in progress provides useful insights into the properties of dusty tori 
around AGN. Fitting the observed optical-to-MIR SEDs to a series of models, 
we showed that in most of the cases the density decreases with distance from 
the centre and can vary with the angle from the equator. Tori properties such 
as the outer radii or the mass of dust are independent of the redshift and only 
depend on the luminosity of the objects. The equatorial hydrogen column density 
distribution is in agreement with what is expected from the Unified Scheme. 
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